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1  IT  ■ 


I.  SEMIANNUAL  REPORT  SUMMARY 

The  projects  being  carried  out  under  this  contract  are  summarized  below. 

More  detailed  discussions  are  given  in  Sections  II-IV  of  this  report. 

(1)  Metal  Vapor-Rare  Gas  Discharges 

The  objective  of  this  project  is  the  evaluation  of  the  potential  of  elec¬ 
trical  discharges  in  high  pressure  metal  vapor-rare  gas  mixtures  for  excitation 
of  high  power,  high  efficiency  lasers  operating  at  near  visible  wavelengths. 

Our  work  has  provided  stimulated  emission  coefficients  for  several  metal  vapor- 
rare  gas  exclmers  at  visible  or  near  visible  wavelengths  and  has  shown  that 
electrical  discharges  with  sufficiently  high  electron  temperatures  (>0.5  eV) 
are  an  efficient  way  to'  produce  these  excimers.  Our  recent  experiments  under 
this  contract  have  shown  that  the  effective  electron  temperature  in  high  power 
discharges  in  Na-Xe  mixtures  with  Na  densities  between  10*"*  and  10*"**  cm  ^  is 
too  low  to  yield  useful  gain  on  the  Na-Xe  excimer  transition. 

During  this  report  period  we  have  completed  work  on  a  detailed  model  of 
the  high  power  discharges  in  Na-Xe  mixtures  which  we  have  examined  experimentally 
under  this  contract.  The  objective  of  this  theoretical  investigation  is  to 
determine  the  dominant  electronic  and  atomic  colllsional  processes  that  control 
the  optical  and  electrical  properties  of  these  discharges.  This  information 
can  then  be  used  to  predict  and  hopefully  to  improve  the  suitability  of  these 
discharges  for  high  power  laser  applications. 

High  power  (~I0®  W/i)  discharges  in  metal-doped  Xe  are  modeled  for 

ic  i  n  20  —3 

typical  metal  atom  densities  of  10  -10  cm  and  Xe  densities  of  ~10  cm 

and  electron  densities  of  lO^-lO1^  cm”'*  as  appropriate  for  proposed  excimer 
lasers.  Na  is  used  as  a  prototype  species,  while  its  properties  are  varied 
to  indicate  some  of  the  changes  that  could  result  from  the  use  of  different 
metals.  The  model  includes  sixteen  excited  levels  of  Na,  three  ions,  and 
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the  excimer  levels  of  NaXe  and  Na2»  Hie  degree  of  ionization  is  determined 
by  collisional  multistep  excitation  and  ionization  of  excited  atoms  vs 
dissociative  recombination  of  electrons  with  Na2+*  Steady-state  conditions  in 
the  positive  column  are  calculated  for  typical  gas  temperatures  of  ~0.06  eV 
and  electron  temperatures  Te  of  0. 3-0. 5  eV.  The  Na  population  distribution  is 
largely  Boltzmann  at  the  electron  temperature  and  the  electron  density  is 
close  to  the  Saha  equilibrium  value  except  at  low  electron  temperatures  and 
very  high  extracted  laser  powers.  Useful  gain  and  extracted  powers  of  ~10 
MW/cm^  are  predicted  for  the  higher  and  Na  densities,  with  the  pulse  width 
limited  to  ~10~^  sec  by  gas  heating.  The  model  indicates  that  a  metal  which 
produces  a  deeply  bound  product  state  via  dissociative  recombination  could 
yield  a  very  efficient,  high  power  laser  or  gain  cell. 

(2)  Electron  Excitation  Rates  in  Metal  Vapors 

The  objective  of  this  new  project  is  to  determine  effective  rate  coeffi¬ 
cients  for  the  electron  excitation  and  deexcltatlon  of  excited  mercury  atoms 
for  use  in  the  modeling  of  the  role  of  these  excited  atoms  in  electrically 
excited,  high  power  lasers,  e.g.,  mercury-halogen  lasers.  In  order  to  minimize 
problems  with  chemical  reactivity,  these  experiments  are  being  carried  out  in 
mercury-rare  gas  mixtures.  The  electrical  and  optical  properties  of  discharges 
in  this  mixture  are  being  investigated  under  another  contract.  This  work  is  a 
prerequisite  to  the  rate  coefficient  determinations  and  is  currently  directed 
toward  (a)  absolute  calibration  of  the  optical  detection  system  so  as  to  allow 
more  accurate  determinations  of  the  densities  of  the  excited  states  of  mercury 
and  (b)  installation  of  an  auxiliary  discharge  so  as  to  provide  initiating 
electrons  and  so  improve  the  stability  of  the  main  discharge.  Once  these 
discharges  are  operating  stably  and  are  characterized,  the  present  investigation 
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wlll  make  use  of  laser  perturbation  techniques  and  fast  transient  analysis  of 
the  emitted  radiation  to  obtain  the  desired  rate  coefficients.  The  fast 
transient  digitizer  and  the  pulsed  dye  laser  for  this  experiment  have  been 
acquired  and  tested.  This  means  that  the  desired  measurements  will  not  be 
made  until  near  the  end  of  this  contract. 

(3)  Electron  Excitation  of  Metastable  Atoms  and  Molecules 

During  this  report  period  we  assembled  and  tested  apparatus  for  measure¬ 
ment  of  electron  excitation  rate  coefficients  for  the  lower  metastable  and 
radiating  states  of  neon.  Electron  excitation  rate  coefficients  for  the  lower 
excited  states  of  the  rare  gases  are  essential  to  the  prediction  and  under¬ 
standing  of  the  efficiency  and  power  output  of  electrical  discharge  excited, 
high  power  rare  gas-halogen  lasers.  In  particular,  recent  work  at  the  Naval 
Research  Laboratory  has  shown  that  the  use  of  Ne  as  the  principal  rare  gas 
results  in  significant  improvements  in  rare  gas-halogen  laser  performance. 

This  result  and  the  convenience  of  neon  wavelengths  led  to  our  choice  of  neon 
for  this  demonstration  experiment.  Resources  will  not  be  available  to  extend 
these  measurements  to  the  other  rare  gases. 

The  technique  being  used  for  the  measurement  of  the  neon  metastable  and 
resonance  state  densities  is  that  of  the  absorption  of  radiation  in  transitions 
from  the  metastable  (or  resonance)  state  to  a  more  highly  excited  state.  In 
order  to  eliminate  the  usual  uncertainties  in  derived  metastable  densities 
caused  by  uncertainties  in  the  spectral  profile  of  the  incident  neon  radiation, 
we  have  replaced  the  usual  neon  discharge  source  with  a  single  mode  dye  laser, 
i.e.,  with  a  source  whose  spectral  width  is  narrow  compared  to  the  absorption 
profile.  The  use  of  a  narrow  laser  beam  also  allows  us  to  use  a  multipass 
absorption  cell  and  so  bring  the  absorption  signal  up  to  levels  well  above  the 
noise  level  of  the  system.  This  arrangement  has  been  tested  and  measurements 
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will  be  made  during  the  next  report  period. 

II.  METAL  VAPOR-RARE  GAS  DISCHARGES 

Drs.  L.  W.  Schumann,  A.  Gallagher  and  A.  V.  Phelps 

We  have  chosen  to  model  the  discharge  excited  Na-Xe  excimer  system  because 
of  the  relatively  large  amount  of  collisional  cross  section  data  available  for 
Na  and  because  of  the  availability  of  experimental  data  for  comparison  with 
our  theory.  It  is  believed  that  many  of  the  properties  of  these  particular 
species  are  characteristic  of  other  mixtures,  as  long  as  the  ionization  energy 
of  the  dopant  is  below  the  first  excited  state  of  the  noble  gas,  and  the  elec¬ 
tron  affinity  of  the  dopant  is  small  enough  that  negative  ion  densities  are 
insignificant.  Much  of  this  work  was  carried  out  earlier  by  R.  Shuker. 

Note  that  for  completeness  this  discussion  includes  some  results  of  the 
measurements  on  TJl-Xe  discharges  carried  out  under  another  contract  (AFWL 
Agreement  78-036) . 

A.  DETAILS  OF  THE  MODEL 

a)  Discharge  conditions 

The  pulsed  discharges  under  consideration  typically  occur  in  Xe  at  densities 

iq  20  15  17 

of  OCe]  *  10  -10  cm  doped  with  Na  at  densities  of  ENa]  *  10  -3*10  cm 

The  electron  densities  (ng)  are  typically  10^-1017  cm  \  the  electric  field 

—18  —17  2 

to  total  gas  density  ratios  E/N  are  in  the  10  -10  V  cm  range,  the  electron 

temperature  Tg  is  0.3-0. 5  eV,  and  the  gas  temperature  T  s  800K,  so  that  kT  s  0.07 
eV.  Here  and  throughout  this  paper  electron  temperatures  are  expressed  in  eV. 

4 

Typical  modeled  discharge  power  densities  are  10-10  MW/Jt.  It  is  assumed  that 
quasi-steady-state  discharges  will  be  terminated  before  excessive  300°C )  gas 

heating  occurs.  The  present  models  thereby  emphasize  steady-state,  constant 
gas  temperature  conditions.  Results  in  Sec.  IIB  are  obtained  by  integration  of 
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nonlinear,  time  dependent  rate  equations,  and  in  Sec.  IID  by  use  of  a  two- 

temperature  approximation.  The  near-uniform  cylindrical  spatial  character  of 

the  discharge  under  consideration  is  described  in  Ref.  1.  Here  we  consider 

only  a  uniformly  excited  region  in  the  positive  column  of  the  discharge. 

The  fractional  ionization  in  the  discharges  is  high  enough  so  that  electron- 

electron  collisions  could  be  expected  to  largely  Maxwellianize  the  electron 

2 

energy  distribution.  Thus,  the  calculations  presented  here  utilize  Maxwellian 
electron  energy  distributions  and  Na  excited  states  up  to  an  effective  quantum 
number  of  15. 

b)  Atomic  states 

An  Na(3P)/Na(3S)  fraction  of  ~0.05,  which  is  equivalent  to  an  ~0.5  eV 
excitation  temperature,  is  desired  for  efficient  laser  action.  (As  another 
example,  T£(7S)/T?„(6Pjy2)  -  0.02  equivalent  to  Tfi  =  0.85  eV  is  desired  for  a 
Ti-Xe  excimer  laser.  )  If  this  excitation  temperature  results  from  a  balance 
between  electron  excitation  and  deexcitation  of  Na(3P),  then  an  electron 
temperature  T£  of  ~0.5  eV  is  necessary.  Electron  collisions  with  Na(3P)  atoms 
will  then  tend  to  produce  a  similar  excitation  temperature  between  Na(3P)  and 
higher  excited  states.  Radiation  will  lower  the  excited  state  populations, 
but  at  the  electron  densities  in  these  discharges  this  is  a  minor  effect.  If 
the  excited  Na  densities  are  characterized  by  an  excitation  temperature  in 
this  0.5  eV  range,  then  the  electron  collisional  excitation  and  ionization 
rates  out  of  excited  levels  increase  with  increasing  excitation  energy  of  the 
level.  Thus,  highly  excited  states  must  be  considered  in  the  ionization 
balance.  We  have  included  Na  states  with  an  effective  principal  quantum 
number  n  up  to  15  in  the  present  model,  and  we  allow  for  the  effect  of  still 
higher  levels  through  the  collisional  rate  coefficients  assigned  to  this 
terminal  n  ■  15  level. 


The  highly  excited  Na  states  are  strongly  mixed  by  electron  collisions 
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so  that  most  of  the  population  resides  in  the  high  angular  momentum  states. 

These  states  are  very  hydrogenic  and  very  nearly  degenerate  at  integral  n 

values.  Consequently  we  have  simplified  the  level  structure  by  grouping  all 

states  above  4P  into  11  hydrogenic  levels,  n  =  4  to  15,  as  shown  in  Figs.  1 

and  2.  This  also  allows  the  use  of  hydrogenic  cross  sections  for  electron 

collisions  with  the  highly  excited  states.  The  rate  coefficients  for  these 

3  4 

collisions  are  obtained  from  the  literature  '  and  by  scaling,  and  are  shown 
in  Fig.  3.  The  lowest  excited  states  of  Xe .  i.e.,  a  single  Xe*,  and  ionization 
(e  +  Xe*  j  2e  +  Xe+)  are  included  in  the  model.  In  fact,  the  Xe+  ions  are 

rapidly  converted  into  Xe£+  ions,  which  undergo  dissociative  recombination  to 

*  +  * 
form  Xe  .  For  the  T  considered  here,  the  densities  of  Xe  and  Xe  are  much 

e 

+  * 

lower  than  Na  and  Na  ,  and  they  do  not  significantly  influence  the  kinetics, 
c)  Molecular  species 

The  lower  energy  levels  and  ionic  forms  of  Na,  and  NaXe  are  shown 

in  Fig.  2.  The  A-X  band  of  the  NaXe  excimer  is  a  broad  continuum  extending 
from  the  590  nm  Na  resonance  line  to  ~720  nm;  the  700  nm  region  of  this  band 
is  the  most  favorable  region  for  excimer  laser  action.  ^  The  B  state  of  NaXe, 
not  shown  in  Fig.  2,  is  a  repulsive  state  correlating  asymptotically  with 
Na(3P)  and  Xe^S^).  The  B-X  band  of  NaXe  extending  from  560-590  nm  is  very 
weak  and  from  thermodynamic  arguments  is  not  expected  to  exhibit  net  gain. * 

The  4SE  state  radiates  a  4SE-XZ  band  ^  which  could  have  net  gain  in  the  present 
discharges  if  sufficiently  high  excitation  temperatures  are  obtained.  Higher 
levels  of  NaXe  are  not  shown,  but  they  are  expected  to  exhibit  a  Rydberg-like 
shape,  approximately  paralleling  the  NaXe+  potential. ^ 

At  the  typical  Na  densities  of  10^-10”^  cm""'*  under  consideration,  a  few 
percent  of  the  Na  is  bound  as  Na,,.  The  associated  A-X  (560-800  nm)  and  B-X 
(450-500  nm)  bands  of  Na^  can  cause  severe  absorption,  although  the  A-X  band 


may  also  be  capable  of  yielding  net  gain  in  a  discharge 
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The  radiative 


1  *4*  3 

lifetimes  of  the  Na~  A  and  B  states  are  about  10  nsec.  The  A  E  ,  an  ,  and 
Z  u  u 

1  *4*  9 

X  Eg  states  of  Na^  are  strongly  bound  so  that  the  Xe  induced  collislonal 

dissociation  rate  coefficient  is  very  small.  The  relative  populations  of  the 

A  and  a  states  are  expected  to  be  maintained  in  approximate  statistical 


equilibrium  and  are  represented  in  the  model  by  an  A  state.  Since 
molecules  are  depleted  by  electron  induced  dissociative  excitation,  the  Naj 
states  are  expected  to  be  present  at  much  less  than  their  equilibrium 


fractions.  In  order  to  evaluate  this  effect,  the  A  and  X  states  of  Naj  are 

included  in  the  calculation,  although  they  are  minor  constituents. 

The  NaXe(A)  state  is  bound  by  only  a  few  kT,  so  that  the  Na(3P)  and 

NaXe(A)  densities  are  maintained  at  their  gas  temperature  equilibrium  ratio, 

20  -3 

typically  10%  as  molecules  for  [Xe]  =  2.7  x  10  cm  ,  by  rapid  association 
and  dissociation,  i.e.,  we  assume  equilibrium  for  the  reaction 


Na(3P)  +  2Xe  <--►  NaXe(A)  +  Xe 

The  dissociation  rate  per  NaXe(A)  molecule  is  about  10*®  sec  *  at  typical 
[Xe]  **  1020  cm'*,  so  that  the  equilibrium  ratio  of  [NaXe(A) ]/[Na(3P) ] [Xe]  *=  K^( 
is  maintained  as  long  as  the  destruction  of  NaXe(A)  as  the  result  of 
deexcitation  by  stimulated  emission  and  of  dissociation  and  deexcitation  by 
electrons  is  slower  than  10*®  sec-*.  This  is  expected  to  hold  at  the  ratio 
ne/[Xe]  <  10  and  the  incident  optical  intensities  considered  here. 
Vibrational  and  rotational  excitation  of  NaXe  A-state  molecules  by  electrons 
is  expected  to  be  slow  compared  to  vibrational  and  rotational  relaxation  by  Xe 
so  that  these  electron  collisions  can  be  ignored.  The  observation  in  Ref.  1 
of  a  nearly  thermal  excimer  band  shape,  corresponding  to  a  vibrational 
temperature  at  the  gas  temperature  T,  also  supports  these  expectations.  In 
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the  present  calculation  we  therefore  include  only  molecular  states  connecting 

to  the  Na(3p)  states,  with  the  implicit  assumption  that  the  equilibrium 

portion  of  this  is  in  the  form  of  NaXe  A-state  molecules.  Thus  all  inelastic 

collision  cross  sections  of  electrons  with  the  molecules  are  effectively 

assumed  to  be  identical  to  those  for  the  parent  Na  state. 

The  Na(3P)  states  radiate  with  a  lifetime  of  16  nsec  while  the  NaXe  A  and 

B  states  have  slightly  longer  and  shorter  radiative  lifetimes,  respectively. 

Therefore,  the  combined  Na(3P)  +  NaXe(B)  +  NaXe(A)  population  radiates  with  an 

effective  lifetime  near  16  nsec.  From  the  normalized  emission  spectrum  of 

this  composite  population,  and  the  absorption  coefficient  derived  therefrom, ^ 

one  finds  about  a  10%  probability  that  a  photon  is  emitted  in  an  optically  thin 

20  —  V 

portion  of  the  spectrum  (for  0.2  cm  escape  distance  and  [Xe]  ~  10  cm  ), 
i.e. ,  the  effective  radiative  lifetime  for  this  combined  population  is  about 

160  nsec.  This  effective  radiative  rate  constant  is  not  very  sensitive  to  [Na] 

2 1  —3 

and  it  scales  linearly  with  [Xe]  for  values  below  10  cm  .  The  stimulated 

“16  2 

emission  coefficient  at  700  run  is  ~10  cm  per  NaXe(A)  molecule,  and  thus  is 

—  17  2  20  — 

~10  cm  per  total  Na(3P)  population  (bound  and  free)  at  [Xe]  =  10  cm 

Some  of  the  transitions  from  higher  Na  states  to  Na(3P)  are  also  slightly 

optically  thick,  but  as  their  radiative  rates  are  much  smaller  than  their 

electron  collisional  rates  this  is  ignored  in  the  calculations. 

The  NaXe+  molecular  ion  is  structurally  an  Na+-Xe  molecule,  whose 

potential  (Fig.  2)  we  take  to  have  ~0.3  eV  binding  and  Rg  ~  3. 2  A  from 

experiments  and  theory. ^  (We  have  no  information  about  T£Xe+. )  The  fact 

that  only  one  molecular  ground  state  results  from  this  interaction  between 

closed  shell  (*Sq)  monomers  has  a  profound  effect  on  its  collisional  and 

optical  properties.  Thus  NaXe+  cannot  absorb  in  the  visible,  and  the  highly 

excited  states  of  Na  -Xe  will  all  nearly  parallel  the  single  Na  -Xe  potential. 


4 


The  latter  property  yields  the  expectation  that  Xe  collisions  will  not 

A  * 

efficiently  quench  highly  excited  Na  atoms  or  bound  Na  -Xe  molecules  and,  as 
discussed  further  below,  that  dissociative  recombination  with  NaXe+  wili  have 
a  small  rate  coefficient.^  The  reactions  Na+  +  2Xe  +  NaXe+  +  Xe  will  be  very 
fast ^  at  high  [Xe] ,  yielding  a  ratio  of  [NaXe+]/[Na+] [Xe]  very  close  to  the 
equilibrium  ratio  at  all  n£  considered  here.  As  discussed  in  Sec.  IIB, 
the  size  of  K*  is  very  important  in  determining  the  effective  rate  of 
electron  recombination  in  the  discharge. 

The  Na2+  molecular  ion  potentials,  from  Ref.  9,  are  shown  in  Fig.  2. 

Under  the  conditions  considered  below  Na2+  is  a  minor  constituent  in  the 
discharge,  so  that  absorption  by  Na2+  (see  Ref.  11)  and  emission  from  excited 
states  of  ttaj  are  expected  to  be  negligible.  However,  e  +  Na2  Na  +  Na  is 
the  dominant  electron  recombination  mechanism  in  our  model.  Na2+  wiH  be 
formed  by  the  exchange  reaction  Na  +  NaXe+  -►  Na2+  +  Xe.  Due  to  the  low 
[Na]/[Xe]  ratio  the  equilibrium  ratio  of  [Na2+]/[NaXe+]  is  typically  on  the 
order  of  1;  rapid  dissociative  recombination  of  Na2+  normally  makes  this  ratio 
much  smaller. 

The  highly  excited  states  of  Na  are  expected  to  have  associated  NaXe 
states  with  binding  energies  about  equal  to  that  of  NaXe  .  If  predissociation 
of  lower  excited  states  of  NaXe  occurs  to  form  excited  Na,  the  loss  of  NaXe* 
will  cause  a  net  downward  flow  of  excitation  In  the  molecular  excited  state 
ladder  and  produce  departures  of  the  Na*  populations  from  Boltzmann  equilibrium 
in  the  same  manner  as  dissociative  recombination  of  electrons  with  NaXe+. 
However,  as  discussed  in  Sec.  IIB,  the  experimental  lower  excited  state 
population  data  seem  to  show  that  departures  from  the  Boltzmann  plot  are 
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d)  Electrons 

At  the  typical  Ionization  fraction  of  interest  to  us,  i.e.,  ne/[Xe]  M.0  \ 

energy  exchange  by  electron-electron  collisions  is  much  more  rapid  than  energy 

* 

loss  by  recoil  of  Xe.  Inelastic  and  superelastic  collisions  with  Na  are 

largely  in  balance  so  that  the  net  energy  transfer  rate  between  electrons  and 
* 

Na  is  much  smaller  than  between  free  electrons.  Therefore  the  electron 

energy  distribution  is  expected  to  be  very  nearly  Maxwellian  at  a  temperature 

Tfi.  The  range  of  validity  of  this  assumption  for  the  higher  excited  state 

12 

densities  is  being  investigated  separately. 

e)  Collision  rate  coefficients 

* 

The  simplified  model  used  considers  the  states  shown  in  Fig.  1,  plus  Xe  , 

Xe+  and  Xe2+.  In  addition  to  Na2(X),  Na^CA),  Na+,  NaXe+,  and  Na2+,  hydrogenic 
levels  of  Na  from  n  =  4  to  15  and  the  3S,  3P,  AS,  3D,  and  4P  levels  of  Na  are 
included.  Electron  collision  rate  coefficients  between  n  **  4  to  15  levels  and 
for  ionization  of  these  levels  are  calculated  from  the  hydrogenic  cross 

3 

sections  of  Johnson,  by  integrating  over  the  assumed  Maxwellian  energy  distri- 

4 

bution,  and  are  shown  in  Fig.  3.  Calculations  by  Moores  et  al  are  used  for 

cross  sections  within  the  five  lower  Na  levels,  and  scaling  and  detailed  balance 

arguments  are  used  for  rate  coefficients  between  the  five  lower  Na  and  the 

n  *  4  to  15  levels.  In  order  to  account  for  the  n  >  15  levels  that  are  neglected, 

we  give  the  n  =  15  level  an  ionization  rate  coefficient  equal  to  about  one  third 

the  rate  for  the  n  =  15  to  n  =  16  step.  This  is  close  to  the  effective  multi- 

13 

step  ionization  rate  coefficient.  The  n  *  14  ionization  rate  coefficient  is 
found  by  interpolation  between  the  n  **  13  and  n  ■  15  values.  Since  the  electron 
energy  distribution  is  Maxwellian,  all  i  -*■  j  rate  coefficients,  k^ ,  are  related 
by  the  thermodynamic  ratio  ■  (gj/g^)expC(E^-E^)/kT^3.  The  thermodynamic 

(Saha)  equilibrium  condition  is  also  used  to  relate  ionization  and  inverse 
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ionization  rate  coefficients.  This  condition  is  usually  written  as 


ne  [Na*) 
8e  8+ 


[Na(n) ] 

8„ 


2xmkT 


3/2 


expJ-EjOO/kT^ 


(1) 


where  Ej(n)  is  the  ionization  potential  of  the  n’ th  level.  The  statistical 

weights  are  gg  =  2  for  the  electrons  and  g+  =  1  for  Na+*  The  gn  vary  from  2 

2 

for  the  3S  level  to  2n  for  the  hydrogenic  levels. 

The  only  heavy  particle  collision  rates  entering  the  rate  equations  are 

g 

those  for  Na2(X)  and  Na2(A')  association  and  dissociation  (as  noted  above 
the  NaXe(A)/Na(3P)  ratio  is  assumed  in  equilibrium),  i.e. , 


Na(3S)  +  Na  +  Xe  t  Na2(X)  +  Xe  +  0.5  eV  , 

and 

Na(3P)  +  Na  +  Xe  X  Na2(A')  +  Xe  +  1  eV  , 
and  for  the  ion-molecule  processes 

Na+  +  2Xe  t  NaXe+  +  Xe  +  0. 3  eV  , 

and 

NaXe+  +  Na  t  Na2+  +  Xe  +  0.  7  eV 

We  assume  that  all  molecules  are  in  an  equilibrated  vibrational  and  rotational 
distribution  at  T,  and  the  ratio  of  forward  to  backward  rate  coefficients  is 
given  by  equilibrium  relations  at  the  gas  temperature  T. 

One  of  the  most  important  processes  controlling  the  behavior  of  these 

discharges,  in  which  most  of  the  ions  are  molecular  and  T  <<  T  ,  is  dissocia- 

e 

tive  recombination  of  electrons  with  the  molecular  ions.  Dissociative 

14 

recombination  is  generally  described  in  terms  of  a  curve  crossing  of  the  disso¬ 
ciative  diatomic  state  (AB*)  with  the  molecular  ion  state  (AB*).  In  homonuclear 
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diatouics  such  crossings  always  occur  because  symmetry  considerations  lead  to 

repulsive  as  well  as  an  attractive  molecular-ion  state;  thus  some  of  the 

Rydberg-like  A2*  states  will  parallel  this  repulsive  state  and  cross  the 

(attractive)  A£+  state.  This  situation  applies  to  Ua^  dissociative 

recombination,  so  that  we  expect  it  to  have  the  typical  large  rate  coefficient 

observed  for  the  heavier  inert  gas  diatomics, 15  e.g. ,  2  x  10”^  (0.026/Te)*^ 

cm  /sec.  For  NaXe  ,  on  the  other  hand,  only  one  attractive  state  occurs  and  the 

Rydberg-like  NaXe*  states  will  not  cross  this  ion  state.  The  rate  coefficient 

for  NaXe+  dissociative  recombination  is  thereby  expected^  to  be  much  smaller 

than  the  Na2+  value.  Including  NaXen+  clusters  does  not  modify  this  conclusion. 

For  the  same  reasons  associative  ionization  to  form  NaXe+  is  assumed  not  to 

occur.  Consequently,  we  have  taken  ct(NaXe+)  =  0  in  most  of  the  calculations. 

4.  4. 

Unless  the  rate  coefficient  for  Na  +  Na  +  Xe  +  Na2  +  Xe  were  usually  large, 

e.g.,  10  cm  /sec  ,  the  Na2  is  formed  predominantly  from  NaXe  through 

the  NaXe+  +  Na  -►  Na2+  +  Xe  ion  exchange  process.  Due  to  the  low  Na  density 

this  process  is  usually  the  rate  limiting  step  that  determines  the  electron- 

ion  loss  rate  in  the  discharge.  For  this  ion  exchange  rate  coefficient  kg  we 

use  a  Langevin  rate  coefficient,  and  believe  that  this  should  be  fairly 

accurate.  However,  the  final  Na*  product  state  of  the  Na2+  dissociative 

recombination  must  be  estimated  by  fitting  the  model  to  the  experiments. 

The  curve-crossing  model  predicts  that  the  final  state  products  of  Na2+ 

dissociative  recombination  (^2"*"  +  e  +  Na(nR)  +  Na)  will  be  an  atomic  state 

bound  by  about  2.3  eV,  the  vertical  dissociation  energy  of  Na2+.  This  would 

apply  to  Na2+  in  the  ground  vibrational  state,  whereas  the  same  model  would 

predict  that  an  excited  vibrational  state  would  dissociate  to  less  strongly 
* 

bound  Na  states.  Thus,  we  expect  a  distribution  of  final  states,  labeled  by 
principal  quantum  numbers  nR,  with  the  mean  values  of  nR  increasing  with  Na2+ 
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vibrational  temperature.  For  convenience  we  often  use  a  single  n_,  representing 

R 

the  effective  value  of  this  distribution. 

Destruction  of  excited  Na  in  collisions  with  Na,  by  associative  ionization 
to  form  Na2+  +  e,  occurs  with  a  typical  gas  kinetic  rate  coefficient.***  Due  to 
the  low  Na  density  it  is  negligible  compared  to  multistep  ionization  by 
electrons.  Similarly,  ion  pair  formation  (Na*  +  Na  ♦  Na+  +  Na  )  is  negligible. 

We  do  not  expect  excitation  and  deexcitation  of  the  highly  excited  Na 
states  in  collisions  with  Xe  atoms  to  compete  with  electron  excitation  and 
deexcitation,  since  these  Rydberg-like  NaXe  potentials  should  be  approximately 
parallel  to  the  single  NaXe+  potential  and  do  not  provide  the  necessary  curve 
crossings.  Even  in  the  case  of  the  homonuclear  system  of  H*  collisions  with  H 
atoms,  theory*^  shows  the  electron  collisions  dominate  for  [e]/[H]  >  We 

have  tested  for  the  influence  of  a  n  =  5+4  quenching  rate  coefficient  of 

3  x  10"10  cm'Vsec-*,  the  maximum  rate  that  appears  feasible  in  the  absence  of 

18 

long-range  curve  crossings,  and  found  it  negligible  under  typical  discharge 
conditions.  Collisional  deexcitation  of  higher  states  is  expected  to  cause 
still  smaller  population  changes,  so  no  deexcitation  by  Xe  is  included  in  the 
calculations  reported  here. 

The  model  described  above  is  the  basis  for  formulating  a  set  of 

differential  equations  which  describes  the  various  processes  occurring  in 

the  discharge  and  the  time  evolution  of  the  species'  densities  and  level 

populations.  The  rate  equations  are  stiff  differential  equations  and  are 

19 

solved  numerically  using  a  Gear  package. 

B.  RESULTS 

The  efficiency  and  pulse  energy  of  lasers  based  on  the  discharges  dis¬ 
cussed  in  this  report  are  greater  if  they  can  operate  during  the  relatively 
long  steady-state  portion  of  the  discharge.  This  is  due  to  the  energy  initially 
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required  to  Ionize  a  major  fraction  of  the  Na  atoms.  The  data  In  Refs.  1  and 
20,  to  which  we  will  compare,  also  represent  steady-state  conditions.  Conse¬ 
quently  we  will  concentrate  here  on  the  calculated  steady-state  densities  of 
the  17  Na  states  in  the  model,  and  of  the  ion  and  electron  densities.  The 
calculations  were  made  for  fixed  electron  temperatures.  As  expected,  there  is 

an  exponential  buildup  of  ng  from  the  assumed  starting  electron  density  of 
10  -3 

10  cm  ,  an  increasing  rate  of  ionization  as  ne  increases,  and  once  ng  exceeds 
14  -3 

about  10  cm  ,  convergence  to  the  steady-state  conditions  within  a  few  nano¬ 
seconds.  In  cases  where  the  neutral  Na  density  is  severely  depleted  at  steady 
state,  the  excited-state  densities  are  larger  during  the  rapid  ionization 
transient  than  for  steady-state  conditions. 

The  excited  state  and  ion  densities  in  the  steady-state  positive  column 
depend  on  T^,  nR,  and  the  effective  dissociative  recombination  rate  per  Na+  ion. 

The  latter  rate  depends  on  the  CNaXe+3/[Na+3  ratio,  which  is  assumed  in  equilib- 

■f  +  + 

rium  at  the  gas  temperature.  Since  Na£  is  formed  by  NaXe  +  Na  -*•  Na^  +  Xe, 

with  rate  coefficient  k  ,  the  effective  Na„+  dissociative  recombination  rate 

a  £ 

per  Na+  is  proportional  to  kfl  CNaXe+3-  The  actual  ratio  CNaXe+3/CNa+3  =  K*  £Xe3 

19 

is  <0.5  at  the  CXeH  <  4.5  *  10  considered  in  Ref.  1,  and  as  a  result  the 
effective  dissociative  recombination  rate  is  too  small  to  significantly  perturb 
the  CNa  3  and  CNa+3  densities,  i.e.,  these  densities  relative  to  tNa3  are  very 
close  to  the  Saha  ratios  at  Tg,  in  agreement  with  the  Na-Xe  discharge  data  in 
Ref.  1.  The  laser  gain  which  results  from  this  limit  is  considered  in  Sec.  IID. 
In  contrast  to  the  Na-Xe  case,  the  Tf.-Xe  discharge  data  of  Ref.  20  does 


not  fit  the  Saha  relation  for  the  higher  n.  This  implies  that  more  rapid 

dissociative  recombination  and  thus  larger  values  of  K+  [Xe]  and/or  a  much 

eq 

lower  energy  for  the  acomic  lervel  resulting  from  dissociative  recombination. 


In  order  to  investigate  this  type  of  behavior  we  have  used  all  other  cross 
sections  as  appropriate  to  Na  and  Xe,  but  increased  K.^  by  a  factor  of  10-20, 


>  i 
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corresponding  to  0.45  eV  rather  than  0.30  eV  of  molecular  ion  binding,  i.e. , 
[NaXe+l  =  n  .  We  have  then  systematically  varied  T  ,  nD,  and  k  .  We  will 
illustrate  the  effect  of  these  parameters,  and  of  including  NaXe+  dissociative 
recombination,  in  the  exanples  shown  in  Figs.  4  and  5.  The  effects  of 
variations  in  the  stimulated  emission  rate  caused  by  external  radiation  are 
shown  in  Fig.  6  for  the  case  of  nR  «  3P.  The  calculated  values  of  extracted 
power  densities  and  laser  amplifier  efficiency  for  the  nR  «  6  (Saha)  case  and 
for  nR  =  3P  are  shown  in  Fig.  7.  Detailed  discussions  of  these  results  are 
presented  in  the  following  subsections. 

a)  Excited  state  densities 

In  Figs.  4-6  we  show  results  of  our  calculations  in  the  form  of  plots 

of  the  log  of  [Na(i)/g^]  versus  excitation  energy.  We  use  2a^  +  1  for  the 

statistical  weight,  neglecting  the  alkali  spin,  since  the  actual  ground  state 

Na  density  then  appears  on  the  abscissa.  An  excited  state  distribution  at  the 

electron  temperature  (Saha  condition)  would  appear  as  a  straight  line  in  these 

figures.  In  all  the  cases  in  Figs.  4  and  5  the  electron  densities  are  high 

enough  so  that  Na  radiative  rates  are  much  slower  than  electron  collisional 

rates  and  the  [Na(i)/g^]  for  the  lower  energy  states  do  nearly  fit  a  straight 

line  at  the  electron  temperature.  However,  above  the  terminal  Na  state  of 

recombination,  nR,  the  population  distribution  frequently  falls  below  such  a 

straight  line.  In  this  region  of  excited  state  energies,  electron  collisions 

are  exciting  Na*  up  a  ladder  of  increasing  n  and  then  ionizing  the  higher  n 

levels.  The  resulting  Na+  forms  NaXe+  and  Na2+,  which  then  dissociatively 

recombine  to  Na(nR)  +  Na(3S)  to  complete  the  cycle  (see  Fig.  1).  Thus, 

competition  between  excitation  and  ionization  by  electrons  and  recombination 

on  Na2+  results  in  population  ratios  [Na(n+1 ) }gn/ (Na(n) ]gn+j  for  n  >  nR  whose 

average  behavior  is  characterized  by  an  excitation  temperature  between  T£  and 

13 

T.  Using  predicted  scaling  based  on  the  Bates'  approximation,  we  find  that  at 
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n  ^/[Na][Xe]  >  10  ^  (for  the  n  =  6,  T  ■  0.3-0. 5  eV  case)  this  excitation 

k  e 

o 

temperature  approaches  Te;  at  much  lower  ng  /[Na][Xe]  it  should  approach  T. 

In  a  gas  discharge  near  equilibrium,  such  as  a  high  pressure  arc,  with  T  ■  T  , 

inverse  ionization  of  Na+  and  dissociative  recombination  of  Na2+  would  be 

balanced  by  the  processes  of  Na  excited-state  ionization  and  Na*  +  Na 

associative  ionization.  However,  in  the  present  discharges  in  cold  gas  the 

Na2+  density  builds  up  so  that  additional  dissociative  recombination  occurs 

and  drives  the  net  ionization  below  the  Saha  value  appropriate  to  Te* 

In  Fig.  4  we  show  the  effect  of  varying  k  and  n„  at  a  fixed  T  =  0.38 

eV.  The  values  [Najp  =  6.3  x  10*  cm"*  and  [Xe]  =  4.5  x  10*^  cm-'*  are  used 

for  the  four  lower  curves  in  Fig.  4;  these  values  are  representative  of  the 

densities  used  in  the  experiments  of  Refs.  1  and  20.  In  the  lowest  (+)  curve, 

-9  3  -1 

nR  =  6  and  kg  «  1  x  10  cm  s  ;  the  nR  value  is  indicated  by  an  arrow.  Note 

that  the  densities  follow  a  Tg  =  0.38  eV  straight  line  for  n  <  6  and  dip 

slightly  below  this  for  n  >  nR.  In  the  next  higher  (•)  case,  and  kg  is 

increased  to  2  x  10  ^  cm*  s-*.  The  dip  below  a  Tg  =  0.38  eV  line  for  n  >  6  is 

more  pronounced  than  in  the  first  case  because  of  the  faster  rate  of  Na2+ 

formation  and  dissociative  recombination.  The  effect  of  changing  nR  can  be 

seen  as  follows.  The  dashed  line  results  from  distributing  the  dissociative 

recombination  equally  between  nR  **  3D,  4P,  4,  5,  and  6  states,  with  kfl  *  1  x 
-9  3-1 

10  cm  s  .  Note  that  the  Increased  difficulty  of  laddering  upward  from 
these  levels  has  caused  a  major  drop  in  ng  and  the  density  of  high-n  states. 

In  the  (0)  case  kfl  =  1  x  10  cm  s  and  nR  «*  4S  is  still  more  strongly 
bound,  causing  ne  and  the  densities  of  states  with  high  n  to  drop  still 
further. 

To  indicate  density  dependences,  the  top  (x)  curve  in  Fig.  4  again  has 
nR  -  6  and  k  ■  1  x  10  ;  cm  b  ,  as  for  the  lowest  (+)  curve,  but  both  (Na]Q 


and  [Xe]  have  been  increased  by  about  a  factor  of  six.  These  two  ■  6  cases 
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dip  below  the  Saha  (straight)  line  by  about  the  sane  amount,  indicating  that 
there  the  density  dependence  is  relatively  minor  (at  constant  [ Na ] q/ I Xe 1 ) • 

This  is  due  to  a  balance  between  the  effective  recombination  rate  (proportional 
to  [Xc ] )  and  the  effective  ionization  rate  (roughly  proportional  to  ne  which 
varies  as  [Nalp). 

The  Naj+  dissociative  recombination  rate  coefficient  kRR  was  varied 
slightly  in  the  different  cases  in  Fig.  A,  but  the  consequences  are  minor 
since  the  NaXe+  +  Na£+  ion  exchange  is  the  primary  recombination  rate  limiting 
step.  The  value  kDR  =  2  x  10-7  (0.026/Tfi)1//2  was  used  for  the  cross  and  plus 
cases  in  Fig.  A,  kQR  -  1  x  10“  7  (0. 026/Te> 1/,Z  for  the  open  circle  case,  kDR  * 

5  x  10~8  (0.026/T^J1/2  for  the  solid  circle  case,  and  kDR  =  2  x  10-8 

(0.026/Te)^2  to  each  state  from  3D  to  6  for  the  dashed-line  case. 

1 6  —3 

In  Fig.  5  we  show  Tfi  dependence,  for  [Na]Q  =  3.5  x  10  cm  and  [Xe]  ■ 

2.7  x  102®  cm as  would  be  appropriate  for  excimer  lasers.  We  have  used 
nR  =  3P  examples  to  illustrate  departures  from  the  Boltzmann  equilibrium  for 
higher  Te,  i.e.,  for  nR  J  3D  the  densities  would  satisfy  a  Boltzmann 
equilibrium  at  all  n  in  the  T£  =  0.52  eV  case.  In  the  nR  «=  3P,  Tg  ■*  0.52  eV 

case  (open  circles  in  Fig.  5)  the  Na  densities  fall  about  a  decade  below  the 

Saha  limit  at  high  n,  while  in  the  T  *  0. A3  eV  case  they  drop  about  5  decades 

below  this  limit.  This  abrupt  change  from  near  Saha  equilibrium  to  severely 

depleted  [Na*]  occurs  for  any  nR;  for  higher  nR  the  critical  Tfi  is  lower.  The 
abruptness  of  this  change  is  due  to  strong  positive  feedback:  an  increase  in 
the  laddering-ionization  rate  coefficient  with  Tg  raises  nfi,  which  further 
enhances  the  ionization  rate.  In  the  final  example  in  Fig.  A  (7)  we  have 
added  rapid  NaXe+  dissociative  recombination  to  show  the  major  effect  this 
has.  Since  in  this  case  the  recombination  does  not  proceed  via  the  relatively 

J.  X 

slow  NaXe  +  Na  +  Na£  +  Xe  process  it  very  effectively  depresses  ng  and  the 
[Ha*]  at  high  n  values. 


The  relative  Na  densities  for  n  -  12-15  in  our  model  are  normally  (for 
ne  >  10  cm  )  determined  by  Te,  since  the  n  +  n  +  1  collisional  mixing  rates 
increase  rapidly  with  increasing  n,  and  greatly  exceed  the  dissociative 
recombination  flux  for  n  >  10.  Thus  all  states  above  n  =  12  closely  follow  a 
Te  slope.  From  the  Saha  condition,  which  must  hold  for  the  very  high  n,  we 
have  from  Eq.  (1) 


r  *  ,  *, 

iNa  /g  ] 

00 


-3/2 

(2irmkTe/h^)  [Na+  ]n^  /2g+ 
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where  [Na*/g*]oo  is  the  n  +  «>  limit  of  [Na(n))/gn  and  g+  =  1  is  the  Na+ 
statistical  weight.  Thus,  the  calculated  ng  and  [Na+]  should  yield  a 
[Na*/g*]tt>  from  Eq.  (2)  that  is  consistent  with  the  n  =  12-15  values  of 

lit 

[Na(n)]/gn.  Using  Eq.  (2)  we  have  calculated  the  [Na  /g  points  at  the  Na 
ionization  limit  in  Figs.  4  and  5,  and  it  can  be  seen  that  these  fit  an 
extrapolation  of  the  n  *  10-15  densities. 

One  can  also  use  the  experimental  ng  from  Ref.  1  and  a  calculated 
[Na  ]/ne  value  to  obtain  an  experimental  value  of  [Na  /g  from  Eq.  (1). 

For  the  [Na]g  =  6.3  x  10*"*  cm  "*  and  [Xe]  •=  4.5  x  10*  cm  ®  case  the  ratio 

X  X  -f 

[NaXe  ]/[Na  1  is  expected  to  be  near  the  equilibrium  value  of  Kg q [ Xe J  =  0.4 

and  [Na2+]  «  ne,  so  that  [Na+]  +  [NaXe+]  +  [Na2+]  “  ne  yields  [Na+]  =  ne/1.4. 

Using  this  relation  in  Eq.  (2)  with  the  experimental  n£  yields  an  experimental 

point  at  the  Na  ionization  limit.  (The  n£  are  calculated  in  Refs.  1  and  20  from 

21 

measured  current  densities  and  E/N  values  and  calculated  electron  mobilities.  ) 
For  the  data  in  Fig.  6  of  Ref.  1  ne  =  2  x  10*"*  cm-"*  and  3.5  x  10*^  cm"* 
respectively  for  the  0.38  and  0.33  eV  cases.  The  resulting  ionization-limit 
points,  at  5.14  eV,  are  ~4  x  10^  cm"*  for  the  Tg  *  0.38  eV  case  and  ~2  x  10® 
cm  for  the  0.33  eV  case,  i.e. ,  the  0.38  eV  point  falls  slightly  below  the 


extrapolation  of  the  0.38  eV  Saha  line  in  the  figure  of  Ref.  1,  while  the 
0.33  eV  point  falls  about  a  factor  of  4  below  the  0.33  eV  Saha  line. 
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In  the  T£  case  we  do  not  have  even  an  estimate  of  K  for  TiXc  ,  although 

*  A 

the  rapid  decrease  in  [T£  / g  ]  above  4.5  eV  in  Fig.  8  of  Ref.  20  suggests  that 
K*  is  much  larger  than  for  NaXe+  and  the  effective  energies  of  the  T£*(nR) 
states  are  near  4.5  eV.  Since  ng  is  an  upper  limit  for  [Tf.+],  one  can  still 
obtain  an  upper  limit  to  (T£  / g  ]  using  [T£  )  «  n  and  the  experimental  n  in 
Eq.  (2).  This  has  been  done  in  Fig.  8  of  Ref.  20,  where  it  can  be  seen  that 
it  confirms  the  existence  of  a  major  falloff  of  the  high-n  densities  relative 
to  the  Saha  condition.  These  T£  data  bear  considerable  resemblance  to  some  of 
the  calculated  curves  in  Figs.  4  and  5. 


b)  Stimulated  emission  and  efficiency 

For  metal-excimer  laser  operation,  one  would  like  to  obtain  a  large 

population  in  the  first  excited  state,  corresponding  to  a  high  excitation 

temperature,  i.e.,  T^  >  0.5  eV  for  the  Na  case.  If  the  higher  levels  are  also 

at  T  ,  however,  most  of  the  neutral  metal  atoms  will  be  ionized.  Ore  solution 
e 

to  this  problem  is  to  increase  [Na]Q;  then  the  Saha  relation  [Eq.  (2)]  forces 

the  ionization  balance  toward  a  smaller  fractional  ionization.  The  medium 

gain  in  the  Saha  limit  is  considered  in  Sec.  IID.  Another  possible  solution 

occurs  if  the  high-n  levels  are  depressed,  relative  to  Boltzmann  equilibrium 

at  Te,  by  dissociative  recombination,  as  in  some  cases  in  Figs.  4  and  5.  In 

order  to  explore  this  latter  possibility  in  detail  we  have  included  stimulated 

emission  in  the  model.  We  continue  to  use  a  K+  of  ~10  times  the  value  for 

eq 

-f 

NaXe  in  order  to  represent  systems  in  which  dissociative  recombination  can 
drive  the  high  Na*  and  n£  below  the  Boltzmann  equilibrium  values.  The 
expected  n^  are  of  course  unknown,  but  states  bound  by  up  to  2.5  eV  are 
expected  in  the  Na  case  and  the  following  results  are  not  very  sensitive  to 
the  choice  of  n^  ■  3P  (bound  by  3  eV)  and  4S  (bound  by  2  eV).  Thus  we  use  n^ 


3P  in  the  following  illustrative  example 
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In  the  two-temperature  model  with  weakly  bound  exc liners,  the  ratio  of 

excimcr  band  stimulated  emission  cross  section  o*.  "  g  /[Na(3P)J  to  absorption 
cross  section  q&  ■  ky/[Na(3S)]  is  given  by  (1/3)  exp  (h(vQ-v)/kT] ,  where  Vq  is 
the  atomic  transition  frequency.  Ihe  open  circle  case  of  Fig.  6  involves  high 
laser  powers,  where  the  absorption  and  stimulated  emission  rates  are  almost 
equal.  For  this  condition  (Ha(3S)]/[Na(3P)]  is  driven  to  the  above  ratio, 
which  corresponds  to  an  Na(3P)  excitation  temperature  of  T  Vq/(vq-v).  For 
700  nm  radiation  and  T  =  750  K  the  quantity  T  Vq/(vq-v)  is  ~0.40  eV.  We  use 
Tg  =  0.52  eV  in  the  examples  to  be  given,  so  the  laser-induced  transition 
between  3S  and  3P  drives  the  [Na(3P) J/(Na(3S)]  population  temperature  from 
0.52  eV  toward  0.4  eV,  thereby  extracting  power  from  the  medium.  This 
transition  in  the  population  ratio  occurs  when  the  stimulated  emission  rate 
equals  or  exceeds  the  electron  quenching  rate  for  Na(3P)  atoms.  One  can  thus 
estimate  the  maximum  power  available  for  extraction  as  about  equal  to  that 

going  into  3P  quenching  by  electrons;  the  following  calculations  confirm  this 
expected  behavior. 

We  have  calculated  excited  state  densities  and  excimer  band  stimulated 

and  absorbed  power  densities  as  a  function  of  incident  flux  I  for  a  700  nm 

on  —  ^ 

wavelength,  Te  =  0.52  eV,  and  [Xe]  =  2.  7  x  10  cm  using  gain  and 
absorption  coefficients  gy/ [Xe] [Na(3P)]  =  5  x  10-38  cm5  and  k^/ [Xe] [Na( 3S) ]  = 
7.5  x  10~40  cm5  from  Ref.  5  for  T  -  750  K.  For  [Na]Q  =  2.7  x  1017  cm-3  the 
resulting  Na*  population  distributions  are  given  in  Fig.  6  while  the  output 
powers  and  efficiencies  are  given  as  a  function  of  incident  radiation  power  by 
the  solid  curves  in  Fig.  7.  (This  incident  power  would  be  inter-cavity  power 
in  a  self-oscillator.)  The  sharp  peaks  in  the  stimulated  power  and  the  abrupt 
increases  in  efficiency  in  the  n^  «  3P  cases  in  Fig.  7  are  due  to  a  sudden 
drop  in  nfi  by  about  an  order  of  magnitude,  as  seen  in  Fig.  6.  This  occurs 
because  the  [Na(3P)]/(Na(3S)3  ratio  is  lowered  by  the  laser  power.  Since 
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ionization  occurs  by  multistep  excitation  from  Na(3P),  ne  is  therefore  forced 
to  decrease,  and  this  also  makes  the  multistep  ionization  less  rapid,  which 
further  decreases  ng.  At  fixed  Tg  the  discharge  can  thus  undergo  a  type  of 
phase  transition,  from  highly  ionized  Na  and  large  ng  to  slightly  ionized  Na 
and  low  ng.  In  the  [Na]g  =  3.5  *  10  cm  '  case,  also  shown  in  Fig.  7,  this 
transition  to  low  ng  occurs  at  a  lower  incident  radiation  field.  Tlius  far, 
our  models  have  examined  this  effect  only  at  a  fixed  Te  corresponding 
approximately  to  a  fixed  discharge  voltage.  If  the  power  supply  had  an 
internal  impedance  comparable  with  or  greater  than  the  discharge  impedance  the 
discharge  voltage  would  rise,  causing  an  increase  in  Tg  and  tending  to  maintain 
n£  near  a  constant  value.  The  net  gain,  extracted  power,  and  efficiency  should 
then  be  larger  than  indicated  in  Fig.  7.  The  maximum  net  gain  coefficient  in 
Fig.  7  is  gy  -  ky  =  0.03  cm  *  with  gv  =  0.05  cm-1  for  [Na]g  =  2.7  x  10*  ^  CD1  \ 

It  is  important  to  note  that  each  excimer  emission  event  transfers  a  fraction 
(vq~v)/v  -  0.19  of  the  photon  energy  into  gas  heating,  whereas  an  absorption 
event  causes  an  equal  amount  of  gas  cooling.  Therefore  the  efficiency  is  not 
severely  limited  by  the  net  stimulated  emission  in  the  excimer  band. 

For  comparison,  we  have  also  indicated  in  Fig.  7  (short  and  long-dashed 
curves)  the  stimulated  power  that  would  be  extracted  if  n^  =  6,  still  with  the 
larger  .  In  this  case  n£  and  all  [Na*]  are  very  near  the  Saha  values. 

These  stimulated  powers  rise  smoothly  with  increasing  incident  flux,  as  ng  is 
nearly  constant.  At  low  fluxes  net  gains  are  approximately  0.015  cm-1  and 
0.006  cm  *  for  [Na]g  values  of  3.5  n  10*^  and  2. 7  x  10*^  cm-^.  At  high 
incident  fluxes  the  calculated  discharge  efficiency  approaches  about  60%, 
with  the  energy  lost  primarily  by  thermalization  of  the  X  and  A  states  of 
NaXe  following  stimulated  emission  and  by  electron  excitation  of  molecular 


ions. 
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The  maximum  efficiency  for  extracted  power  in  the  nR  =  3P  cases  of  Fig.  7 

is  also  ~60%,  with  19%  going  into  gas  heating  during  rethermalization  of  the 

NaXe  molecules  and  with  the  remaining  ~20%  going  primarily  into  dissociative 

recombination  of  electrons  and  Na2  ions.  Ihis  large  efficiency  is  the  result 

of  the  large  rates  of  Na(3P)  deexcitation  collisions  and  of  the  consequent 

large  rate  of  saturated  stimulated  emission  compared  to  the  effective  rate  of 

2 

ionization  for  Na(3P)  atoms.  The  maximum  extracted  power  at  ~350  MW/cm 

3  3 

incident  laser  flux  is  ~10  MW/cm  or  10  J/£  in  an  0.1  psec  pulse.  (The  pulse 

length  is  limited  to  ~0. 1  psec  by  gas  heating.)  We  must  reemphasize  that  this 
model  applies  to  a  hypothetical  excimer  with  K  >  10  cm  and  in  which  the 
equivalent  of  e  +  Na2+  -*•  Na(3P)  +  Na(3S)  occurs.  Furthermore,  as  discussed  in 
Sec.  IIC,  the  model  does  not  contain  an  unknown  process  that  is  responsible 
for  approximately  75%  of  the  discharge  energy  reported  in  the  measurements  of 
Refs.  1  and  20. 

C.  DISCHARGE  ENERGY  BALANCE  AND  V-I  CURVES 

Once  the  kinetics  have  been  established  (see  Sec.  IIB) ,  one  expects  to  be 
able  to  write  down  the  energy  balance  for  the  metal-rare  gas  discharge.  In 
principle  the  energy  balance  can  be  used  to  derive  volt-ampere  curves  for  the 
discharge.  However,  we  will  show  that  the  energy  loss  processes  considered  in 
our  model  do  not  add  up  to  the  experimental  energy  input  so  that  there  must  be 
energy  loss  mechanisms  that  we  have  not  considered.  Only  the  dominant  energy 
loss  processes  are  discussed  here. 

The  electron  energy  balance  equation  relevant  to  the  present  problem  is 
JE  =  [hvRv  +  vuk(Te-T)ne  +  (eR+3kTe/2)aDRne[Na2+] 

+  (eR’+3kTe/2)a’DRne[NaXe+]  +  ^  jk  ±  NaXe+]  . 


(3) 
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Here  the  power  input  to  the  electrons  is  JE,  where  J  is  the  current  density 
and  E  the  electric  field.  The  first  term  on  the  right-hand  side  is  the  net 
radiative  energy  loss,  the  second  term  is  due  to  elastic  recoil  of  Xe  atoms, 
the  third  and  fourth  terms  are  due  to  dissociative  recombination  of  N'a ^  and 
NaXe+,  and  the  last  term  represents  vibrational  and  rotational  excitation  of 
NaXe  .  is  the  rate  (cm  s  )  of  escape  of  photons  of  energy  hv  and  vu  is 
the  frequency  of  energy  exchange  collisions  of  electrons  of  temperature  Te 
with  the  neutral  atoms  of  temperature  T.  Also,  eR  and  are  the  average 
ionization  energies  of  excited  atoms  produced  by  dissociative  recombination  of 

electrons  and  Naj+  and  NaXe+  ions  respectively  with  rate  coefficients  and 
“dr'  respectively.  cij  and  are  the  energy  transferred  and  rate  coefficient 
associated  with  electron  impact  excitation  of  NaXe+  from  vibrational- 
rotational  level  i  to  level  j.  Note  that  the  first,  third  and  fourth  terms  on 
the  right-hand  side  of  Eq.  (3)  replace  the  large  number  of  terms  representing 
the  electron  energy  loss  and  gain  in  excitation  and  deexcitation  collisions 
with  metal  atoms. 

The  radiative  loss  term  in  Eq.  (3)  includes  spontaneous  emission  and 

stimulated  emission  and  absorption  via  the  excimer  band.  Since  we  have  assumed 

that  the  upper  and  lower  molecular  states  of  the  excimer  transition  are 

maintained  in  thermal  equilibrium  with  the  respective  atomic  states,  the  energy 

supplied  by  the  electrons  per  excimer  photon  emitted  is  equal  to  the  excitation 

1,20 

energy  of  the  atomic  resonance  state.  Experiments  and  our  calculations  show 
that  the  contribution  of  spontaneous  radiation  for  wavelengths  between  200  and 
800  nm  is  only  a  few  percent  of  the  discharge  power,  and  we  believe  that 
radiative  losses  are  still  negligible  when  shorter  and  longer  wavelengths  are 
included.  As  noted  in  Sec.  IIB,  for  Na  the  net  gas  heating  is  19Z  of  the 
optical  energy  extracted  at  700  nm. 
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The  rate  coefficient  for  energy  exchange  collisions  in  Eq.  (3)  is  roughly 
equal  to  (2m/M)<omv> [Xe]  for  electron  temperatures  of  interest  in  this  paper, 
i.e.,  the  electrons  lose  energy  by  elastic  recoil  of  the  Xe  atoms.  Here  m  and 
M  are  the  electron  and  Xe  atom  masses  and  <oEv>  [Xe]  is  the  average  frequency 
of  momentum  transfer  collisions  of  electrons  with  Xe  atoms.  Using  the  cross- 
section  data  from  Ref.  21,  we  obtain  a  normalized  energy  exchange  collision 
frequency  vu/[Xe]  =  1.5  x  10~12  (Te)1,6  cm3/sec_1  for  Tg  between  0.3  and  1  eV. 
We  then  find  that  for  the  experimental  ng  =  2  x  10  cm  of  the  0.38  eV  case 
in  Ref.  1  ,  the  energy  lost  by  elastic  recoil  is  about  10%  of  the  experimental 

power  input  of  about  16  kW/cm  . 

Since  the  Ua2+  dissociative  recombination  rate  is  limited  by  the  ion 
exchange  rate  k  ,  the  contribution  of  dissociative  recombination  to  the  energy 

a 

balance  can  be  estimated  as  ka(Na] [HaXe+]AE,  where  the  energy  loss  per 
recombination  AE  »  3kTe/2  +  cR,  i.e., 

«n»nJN0  =  k  [Na(3S)][NaXe+]  .  (A) 

OK  e  L  a 

Using  a  mean  value  of  ~1  eV  for  eR,  [NaXe+]  =  neK^[Xe]/(l  +  K^[Xe])  and 

_q  -i  _  i 

k  =  10  cm  s  yields  recombination  energy  losses  which  are  typically  about 

o  ' 

A%  of  the  experimental  energy  input  for  the  0.38  eV  data  of  Ref.  1.  (In  the 

limit  of  K+  [Xe]  »  1,  where  K+  refers  to  T£Xe+,  the  recombination  energy 
eq  eq 

loss  is  about  20%  of  the  total  energy  input  for  the  T£-Xe  discharge  data  of 
Ref.  20).  If  NaXe+  dissociative  recombination  to  a  high  nR  (e.g. ,  8)  is 
added,  the  energy  loss  via  dissociative  recombination  is  increased  by  a  large 
factor  without  significantly  changing  nfi  or  the  measured  Na*  values.  Thus 
adding  NaXe+  recombination  at  high  n,  typically  with  a  rate  coefficient  <j|jR 
of  about  half  the  Na2+  value,  is  one  way  to  explain  the  measured  powers. 
However,  as  indicated  in  Sec.  HA  theoretical  arguments  predict  agR  of  NaXe+  to 


be  much  smaller  and  this  contribution  to  the  energy  balance  to  be  negligible 
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(Tliis  also  applies  to  T£-Xe,  and  to  any  metal  whose  ionization  energy  Is  less 
than  the  Xe  excitation  potential.) 

Based  on  the  theory  of  Ref.  23,  rotational  and  vibrational  excitation 
of  the  dominant  molecular  ion  NaXe+  appear  to  be  a  small  (~2%)  but  not 
negligible  energy  drain.  Other  energy  loss  processes  considered  but  rejected 
as  too  small  to  be  of  importance  in  Eq.  (3)  were  a)  collisions  between  two 
excited  atoms;  b)  inverse  ionization  and  collisional  deexcitation  down  the 
NaXe  molecular  excited  state  ladder  to  the  NaXe(A2£)  state  followed  by 

dissociation  and  then  by  electron  induced  excitation  and  ionization  up  the 

£ 

atomic  excited  state  ladder;  c)  quenching  of  Na  states  due  to  Xe  collisions 

with  a  rate  coefficient  of  about  10~^  cm^/sec  *;  d)  dissociation  of  molecules 

v  2A 

by  electron  impact;  and  e)  free-free  radiation,  which  has  been  calculated 

and  also  could  be  detected  at  wavelengths  below  1  pm  if  it  were  large  enough 

to  yield  a  significant  energy  loss.  We  therefore  calculate  a  total  energy 

loss  rate  for  the  Te  -  0.38  eV  case  of  Ref.  1  of  less  than  25%  of  the 

measured  power  input.  A  similar  discrepancy  is  found  for  the  Te  =  0.33  eV 

case. 

An  alternative  means  of  examining  the  energy  balance  is  to  derive  an 
expression  for  the  electric  field  versus  current  density,  or  more  commonly, 
the  volt-ampere  curve  for  the  discharge.  We  first  note  that  the  current 
density  is  given  by  the  relation  (N=[Xe)) 

J  -  en  p  N(F,/N)  ,  (5) 

e  e 

where  pfiN  is  the  electron  mobility  normalized  to  unit  density.  When  Eqs.  (A) 
and  (5)  are  substituted  in  Eq.  (3),  with  radiation  neglected,  one  obtains  the 
equation 

ep  N(E/N)2  -  (cp+3kT  /2)k  [Na(3S)]/N  +  (v  /N)k(T  -T) 
c  k  e  a  u  e 

+  {<eR'+3kTc/2)«DR’  +  ^ljkij}(J/N){epeK(E/N)l 


(6) 
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From  Eq.  (6)  one  sees  that  when  the  energy  loss  is  due  only  to  Na£+  dissociative 

recombination  and  Xe  elastic  recoil  [first  two  terms  on  the  right  side  of  Eq.  (6)3, 

E/N  depends  on  the  discharge  current  density  only  through  the  rather  slow  changes 

1  20 

in  electron  temperature  with  J.  Experimental  measurements  ’  show  that  E/N 

increases  much  more  rapidly  with  J  than  predicted  from  the  temperature  increase 

of  these  terms.  When  the  last  term  in  Eq.  (6)  dominates,  the  E/N  varies  approxi- 

mately  with  (J/N)  .  Approximate  agreement  in  magnitude  and  J/N  variation  with 

the  initial  voltage-current  data,  before  measurements  of  cathode  fall,  led  us  to 

22 

the  no  longer  tenable  hypothesis  that  dissociative  recombination  of  electrons 
with  NaXe+  was  the  dominant  energy  loss  mechanism  in  the  experimental  discharges. 
However,  this  now  appears  to  be  very  unlikely  as  are  other  possible  causes  we 
have  investigated,  such  as  rotational  and  vibrational  excitation  of  NaXe+. 

D.  GAIN  COEFFICIENT  IN  THE  SAHA- BOLTZMANN  LIMIT 

The  primary  objective  of  this  section  is  to  explore  the  consequences  of 

i  •fg 

assuming  that  [Na  ],  ng,  [Na],  and  all  [Na  ]  are  in  Saha- Boltzmann  equilibrium 
ratios  at  temperature  T0.  As  in  the  previous  sections  we  assume  that  atomic 
translation  and  molecular  vibration  temperatures  are  at  the  gas  temperature  T, 
so  that  gain  in  the  excimer  band  results  at  a  sufficient  excited  state  density, 
as  in  the  two-temperature  model  of  Ref.  22.  As  suggested  in  Sec.  IIB  ,  a  rapid 
improvement  in  laser  gain  with  increasing  (Na]Q  is  predicted.  In  the  Tfc-Xe 
case,  the  data  indicate  that  less  ionization  actually  occurs,  perhaps  due  to  a 
larger  binding  of  T£Xe+  compared  to  NaXe+.  A  modification  to  this  simple  two- 
temperature  model  for  this  case  will  be  discussed  at  the  end  of  this  section. 

As  was  found  in  the  detailed  model,  we  assume  that  Na  +  2Xe  «-  NaXe  +  Xe 
is  in  equilibrium  at  the  gas  temperature,  l.e. , 

■ig£*£-I  ■  [Xe)K+  (T)  ,  (7) 

[Na+]  eq 

where  we  use  Kgq(T)  **  6  x  10”^  (T/670  exp(0.3  eV/kT)  cm”"*.  Ftom  the 

calculations  of  Sec.  HB  we  find  that  [ Na 2"** ] / 1 Na"*" ]  «  1,  as  expected  due 
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to  slow  formation  and  rapid  depletion  of  Na2  at  ng  >10  cm  .  Thus,  using 
Eq.  (7), 

n  -  [Na+j  +  [NaXe+]  +  (Na,+  ]  5  [ Na+J{ l+[Xe ]K*  (T)}  .  (8) 

Combining  Eqs.  (1)  and  (8),  we  obtain 


n2  S  <{l+[Xe]K^  (T)}[Na(3S)J(2irmkT  /h2)  exp(-E_/kT  )~  ^  ,  (9) 
e  ecj  e  l  e 

If  we  take  into  account  neutral  depletion,  i.e. ,  [Na(3S))  «  [Na]Q-ne~(Na2+]  » 
[Na]Q-ne  then  the  electron  temperature  is  given  approximately  by  the  solution 
of  the  relation 


-Ej/kT 

e 


e 


n2/( [Na] 0~ne) 

(2T!kTe/h2)3/2{l+[Xe]K^(T)} 


(10) 


Assuning  equilibrium  between  the  Na(3P)  and  Na(3S)  densities  at  Te,  we  see 
that  the  [Na(3P)]  density  can  then  be  expressed  as 


[Na(3P) ] 


exp[-E(3P)/kT  ]  =  [exp(-E  /kT  )]2,1  eV/5‘ 14  eV  #  .  (11) 

6  X  C 


When  Eq.  (11)  is  substituted  into  the  expression  for  the  net  gain,  i.e., 
G  ■  a_[Na(3P)]  -  o„(Na(3S)]  one  obtains  the  relation 

©  SL 


[Na]0(l-X) 


2a 


d-x) 


(12) 


where  x  ■  ne/[Na]0,  B  =■  3[Na]g(  (2?tmkTe/h2)3/2  x  (l+tXejKgq)  ].“%  and 
o  *  2.1/5.14  ■  0.41.  Since  B  varies  slowly  with  T  ,  most  of  the  variation  in 
gain  with  electron  temperature  has  been  absorbed  into  the  variation  with 
electron  density.  We  therefore  maximize  the  gain  in  Eq.  (12)  with  respect  to 
electron  density  at  constant  B  to  obtain  the  maximum  gain  values  shown  in  Fig. 
8.  Note  that  nfi  is  a  controllable  experimental  parameter  determined  primarily 


by  the  discharge  current  density.  At  high  radiative  fluxes  the  Na(3P)  and 
Na(3S)  densities  are  no  longer  in  equilibrium  and  one  must  use  a  more  detailed 
model,  such  as  that  of  Sec.  IIB. 

Figure  8  indicates  that  very  large  gain  coefficients  might  be  achieved  at 

higher  Na  densities  than  those  investigated  in  Ref.  1  (10^-10^  cm-'*). 

However,  two  major  problems  with  the  higher  density  regime  are  the  rapid  gas 

heating  and  the  low  discharge  impedance.  Using  the  relations  given  in  Sec.  IIC 

for  the  energy  dissipation  and  the  conditions  given  in  Fig.  8  for  an  optical 

gain  of  0.02  cm~*,  we  calculate  rates  of  rise  of  gas  temperature  of  about  2000 

K/psec.  The  corresponding  discharge  resistivity  is  0.07  ohm-cm.  Because  of 

3 

the  unexplained  excess  energy  losses  in  the  experimental  discharges  discussed 

in  Sec.  IIC  the  actual  rate  of  rise  of  gas  temperature  and  the  discharge 

impedance  could  be  significantly  larger  than  these  values.  Thus,  while  the 

17  —3 

gain  and  excitation  temperature  are  very  favorable  at  ng  »  [Na]  »  10  cm 
the  power  dissipation  would  limit  the  discharge  pulse  length  to  about  100 
nsec.  Of  course,  this  is  just  a  statement  of  a  universal  problem:  Any  device 
with  an  efficiency  e  for  laser  power  extraction  and  an  allowable  gas  heating 
of  AT  will  be  limited  to  pulse  energy  densities  of  CpNAT  and  extracted  energy 
densities  of  eCpNAT,  where  Cp  is  the  specific  heat  and  N  the  density  of  the 
gas. 

The  Tjl  doped-Xe  experiment  in  Ref.  20  indicated  that  [Ti+]  and  the  higher 

[T£*]  are  below  the  Saha  or  equilibrium  values  at  Tfi  (see  Ref.  20,  Fig.  8). 

According  to  our  model  this  indicates  that  K+  [Xe]  »  1  for  TiXe+.  This 

eq 

reduction  of  [T£+]  below  the  Saha  value  by  a  factor  F  <  1  can  be  accounted  for 

in  Eqs.  (1),  (7)-(ll)  by  multiplying  the  right  side  of  Eq.  (1)  by  F.  The 

quantity  (1  +  K+  [Xe])  in  Eqs.  (9),  (10),  and  (12)  is  then  multiplied  by 
eq 

F.  In  Sec.  IIB  we  tested  for  the  effect  of  such  an  increase  in  K+  by 

eq 


-31- 


increasing  that  for  NaXe  by  a  factor  of  10-20.  Some  reduction  in  F(1  +  K' 


>+ 

eq 


tXe))  did  occur,  but  the  consequence  in  Eq.  (12)  was  less  than  502  reduction 
in  the  gain  shown  in  Fig.  8.  This  is  in  part  because  increasing  K*  increased 
the  effective  dissociative  recombination  rate  and  thereby  decreased  F.  Thus 
the  optimum  gain  appears  to  be  almost  independent  of  the  size  of  K*  ,  although 
the  energy  dissipated  in  the  ionization-recombination  loop  is  directly 


proportional  to  K 


eq 


E.  CONCLUSIONS 


We  have  modeled  a  very  high-power  density  electric  discharge  in  a 

nonequilibrated  vapor  consisting  of  low  ionization-potentials  atoms  (Na  or  T£) 

doped  into  a  high  ionization-potential  atomic  gas  (Xe).  It  was  necessary  to 

include  a  large  number  of  excited-state  species,  plus  several  ion  and  molecular 

species.  At  the  lower  Tg  considered  the  principal  mechanism  controlling  the 

discharge  impedance  is  wultistep  ionization  from  a  state  of  principal  quantum 

number  nR,  followed  by  molecular  ion  formation  and  dissociative  recombination 

back  to  nR.  Comparison  to  Tf,  excited-state  densities  in  a  Tg-doped  Xe  discharge 

Indicates  that  nR  values  as  low  as  3  occur  for  this  case.  In  the  Na-Xe 

discharge  case  it  is  not  possible  to  infer  nD  because  the  small  K+  causes 

k  eq 

relatively  small  departures  from  Saha  densities  at  Tg.  In  the  Na  case  the 
measured  discharge  power  typically  exceeds  that  in  the  model  by  a  factor  of 
four.  Electron-Xe  elastic  recoil,  Na2+  dissociative  recombination,  radiation, 
quenching  of  excited  Na  and  NaXe  states,  and  rotational  and  vibrational 
excitation  of  NaXe+  were  considered.  As  yet  we  have  no  explanation  for  this 
extra  power  dissipation. 

At  the  Te  >  0.5  eV  appropriate  for  NaXe  excimer  laser  operation,  the 

densities  of  ng  and  all  Na  states  are  very  nearly  in  Saha  equilibrium  at  Te> 

regardless  of  nR.  Reasonable  gain  coefficients  and  efficiencies  can  be  obtained 

17  -3 

at  high  CNs]q  under  these  circumstances,  e.g.,  for  CNs]q  ■  2.7  *  10A  cm  J  and 
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QXeJ  ■  2.7  «  10  cm  a  maximum  gain  of  0.03  cm  and  delivered  power  of 

3  2 

10  MW/ cm  at  an  incident  radiative  flux  of  350  MW/cm  and  an  electron  density 

17  -3 

of  10  cm  .  A  10  cm  length  of  such  a  medium  could  readily  yield  a  self- 
excited  oscillator,  although  the  power-buildup  e-folding  time  of  'vl  nsec  in 
this  low-gain  medium  must  be  considered.  This  conclusion  does  not  depend  on 
an  understanding  of  all  the  energy-loss  mechanisms  in  the  discharge,  but  the 
efficiency  does.  The  calculated  maximum  efficiencies  of  'V'60%  are  probably  at 
least  a  factor  of  3  too  high,  but  as  we  have  not  identified  the  dominant 
energy-loss  mechanism  or  its  dependence  on  n&  and  T^,  we  cannot  reliably 
estimate  the  efficiency.  The  scaling  of  this  gain  and  power  output  with  UNa^ 
and  other  parameters  allows  predictions  of  operating  conditions  for  maximum 
gain  to  be  made  rather  simply. 

By  considering  the  case  of  a  large  K+  value  and  dissociative  recombination 

eq 

of  Na2  to  produce  Na(3P)  atoms,  we  have  identified  an  interesting  sort  of 
phase  transition  in  which  ng  decreases  and  CNa]  increases  by  about  an  order  of 
magnitude  at  a  critical  laser  power.  In  addition  to  raising  the  discharge 
Impedance,  this  lower  ne  condition  would  probably  yield  higher  efficiency  in 
a  real  case  since  experimental  data  indicate  that  the  unknown  energy-loss 
mechanism  scales  as  a  fairly  high  power  of  ng.  These  results  point  to  the 
desirability  of  finding  a  metal  species  M  with  relatively  strong  MXe+  binding 
for  which  the  dimer  ion  dissociatively  recombines  into  a  strongly  bound  excited 
atom  and  a  ground  state  atom. 


III.  ELECTRON  EXCITATION  RATES  IN  METAL  VAPORS 
Drs.  E.  R.  Mosburg,  Jr.  and  A.  V.  Phelps 

This  new  experiment  is  intended  to  yield  effective  rate  coefficients  for 
electron  collisions  with  highly  excited  atoms  so  as  to  test  models  of  the  high 
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power  density  discharges  in  metal  vapor-rare  gas  for  laser  applications.  We 
propose  to  use  a  tunable  dye  laser  to  modulate  the  density  of  selected  highly 
excited  levels  of  mercury  atoms  in  a  discharge  so  as  to  obtain  effective 
electron  collision  and  ionization  rrft^bo&fWfcSfents  for  these  excited  states. 

As  indicated  in  the  Report  Summary,  work  on  this  portion  of  this  contract 
during  the  current  period  was  limited  to  the  acquisition  of  essential  equipment. 

IV.  ELECTRON  EXCITATION  OF  METASTABLE  ATOMS  AND  MOLECULES 
Drs.  K.  Tachibana  and  A.  V.  Phelps 

During  this  report  period  we  have  designed  and  constructed  an  apparatus 

3 

for  the  measurement  of  electron  excitation  rate  coefficients  for  the  P ^  and 

3 

P^  metastable  and  resonance  states  of  Ne.  A  schematic  of  the  experimental 

apparatus  is  shown  in  Fig.  9.  The  electron  drift  tube  is  the  same  as  that 

used  in  the  measurement  of  excitation  rate  coefficients  for  the  b^E  raetastable 

3  3 

state  of  0^  and  the  A  E  and  C  II  states  of  and  is  described  In  detail  in 
Ref.  25.  Photoelectrons  released  from  the  semitransparent  cathode  are 
periodically  heated  and  caused  to  flow  across  the  drift  region  by  a  square 
wave  voltage  applied  to  the  anode.  The  electrons  excite  the  ambient  Ne  atoms 
to  the  metastable  and  resonance  states.  These  atoms  are  destroyed  by  collisions 
and  lost  to  the  electrodes  by  diffusion  of  Che  metastables  and  by  radiation 
from  the  resonance  atoms. 

The  time  variation  of  the  density  of  neon  metastable  and  resonance  atoms 
is  determined  from  measurements  of  the  absorption  of  narrow  band  radiation 
(200  MHz)  emitted  by  a  single  mode  dye  laser.  The  dye  laser  is  stabilized  to 
the  peak  of  the  neon  absorption  line  using  the  optogalvanic  effect  from  a 
commercial  Ne  hollow  cathode  discharge  lamp.  As  shown  in  Fig.  9,  the  optical 
system  includes  a  multi-pass  (normally  32  passes)  cell  which  gives  an  active 


absorption  path  of  nearly  2  m.  Figure  10  shows  one  of  the  better  absorption 

transients  obtained  thus  far.  In  this  case  the  quasi-steady-state  mctastable 

8  3 

density  corresponds  to  about  10  metastables/cm  and  the  rms  noise  signal 

6  3 

corresponds  to  about  10  metastables/cm  . 

A  second  example  of  the  operation  of  this  system  is  shown  in  Fig.  11, 
where  the  relative  fluorescent  intensity  is  plotted  as  a  function  of  the  laser 
frequency  as  the  laser  is  tuned  through  the  absorption  line  profile.  The  solid 
dashed  curve  shows  the  profile  calculated  assuming  the  expected  Doppler  broadening 
at  room  temperature  and  taking  into  account  the  known  isotope  components.  This 
data  illustrates  the  spectacular  improvement  in  the  accuracy  of  absorption 
measurements  expected  with  this  new  technique.  Instead  of  an  unknown  emission 
profile  from  the  conventional  neon  capillary  discharge  source  and  a  poorly  known 
overlap  with  the  absorption  profile,  we  can  now  use  the  narrow  and  tunable  laser 
source  to  measure  the  fluorescence  profile  and  so  obtain  the  absorption  profile. 

During  the  next  report  period  we  will  use  this  technique  to  measure  metastable 
and  resonance  atom  absorption  transients  over  a  range  of  Ne  densities  and  E/N 
values  so  as  to  obtain  the  data  necessary  for  the  determination  of  the  electron 
excitation  rate  coefficients  for  these  excited  states  of  Ne. 


Figure  1  -  Energy  levels  and  radiative  and  collisional  processes  used  in  the 

model.  Dashed  lines  are  due  to  gas  collision  processes,  solid  lines 
to  electron  collisions,  single  wavy  lines  to  spontaneous  emission, 
and  double  wavy  lines  to  absorption  and  stimulated  emission. 
Ionization  from  and  Inverse  ionization  to  all  Na  levels  are  included 
in  the  model,  although  only  those  to  higher  levels  are  indicated 
here.  Some  vibrational  levels  of  Na2+  are  indicated,  although  Na2+ 
is  treated  as  a  single  state  in  the  calculations.  This  figure  shows 
the  alternate  paths  for  clectron-Nap.4"  dissociative  recombination 
leading  to  Na  atoms  in  the  3P  or  n  ■  6  levels. 


ENERGY  (eV) 


-36- 


gure  2  -  Electronic  energies  of  various  NaXe,  Naj,  and  Na  states. 

The  Naj,  Na2+,  and  Nag-  potentials  are  from  Ref.  9  ,  NaXe+ 
from  Ref.  7  ,  NaXe  and  NaXe(A2n)  from  Ref.  5,  and  NaXe(4E) 
from  Ref.  6. 
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Calculated  Na  excited-state  densities  [Na(i)J  divided  by  g£  =  21^  +_1 
versus  excitation  energy.  The  conditions  are  [Na]o  =  3.5  x  1016  cm  3 
and  [Xe]  =  2.7  *  102°  cm"3  for  the  top  curve,  and  [Na]o  =  6.3  *  1015 
era'3  and  [Xe]  =  A. 5  *  1019  cm"3  for  the  other  cases,  with  Tg  =  0.38  eV 
in  all  cases.  The  states  nR  produced  by  Na2+  dissociative  recombination 
are  indicated  by  arrows.  For  the  top  and  bottom  (x  and  +)  curves  nR 
*  6,  and  ka  =  10“9  cm3  s'"1,  resulting  in  ne  =  l.A  x  1016  and  3.5  x 
1013  cm-3  respectively.  The  other  (•)  np  =  6  case  uses  ka  =  2  x  10-9 
cm3  s"1,  with  the  larger  ka  causing  a  decrease  in  ionization  (ne  = 

2.6  x  1015  cm-3)  relative  to  the  (+)  case.  In  the  (0)  case  n^  =  4S 
and  kfl  =  10-9  cm3  s'"1,  resulting  ng  =  1.1  x  1014  cm~3.  The  dashed 
line  case  uses  the  same  ka  combined  with  Na2+  dissociative  recombination 
distributed  equally  to  the  5  states  between  3D  and  6,  resulting  in 
n.  =  7. A  x  1014  cm-3. 
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Figure  5  -  [Na(i)]/g^  for  [Na]g  =  3.3  *  1016  cm-3,  [Xe]  =  2.7  *  1020  cm-3, 
and  as  a  function  of  the  indicated  electron  temperatures. 

Naj+  dissociative  recombination  to  nR  =  3P,  as  indicated  by  the 
arrows,  is  included  with  knR  =  2  *  10-7  (0.026/T-,) l/2  and 
ka  *  10-3  cmJ  s  .  NaXe"T’  dissociative  recombination  to  nR  =  3P, 
with  the  same  kpR,  is  a] included  in  the  upper  (7)  0.52  eV  case. 
The  resulting  ne  are,  from  top  to  bottom,  1.1  *  1015,  3.5  *  1014, 
and  3.2  x  10* 6  cm-3.  The  points  at  the  ionization  limit  represent 
the  electron  and  ion  density  through  Eq.  (2). 
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Figure  6  -  Same  quantities  as  in  Fig.  4  as  a  function  of  700  nm  radiation 
intensity  I,  with  [Na]g  =  2,7  *  1017  cm-3,  [Xe]  =  2.7  x  102  cm 
Te  -  0.52  eV,  and  n«  -  3P.  The  closed  circles  are  for  I  =  0,  the 
x  for  I  =  320  MW/cm?  and  the  open  circles  for  I  =  560  MW/cm2;  the 
corresponding  values  of  ng  are  2.0  x  1017,  1.65  x  1017,  and 
1.4  x  1016  cm3,  respectively. 


;ure  7  -  Stimulated  radiative  power  extracted  from  the  discharge  (lower 
four  curves)  and  laser  amplifier  efficiency  (upper  three  curves) 
as  a  function  of  incident  optical  power  at  700  nm.  The  solid 
curves  are  For  the  conditions  of  Fig.  6.  The  chain  (  **) 

curves  are  for  [Na)o  *>  3.5  *  10 16  cm-3  and  ng  •*  3P.  The  long 
and  short  dashed  curves  are  for  n^  =  6  and  for  [Na]o  ■  3.5  *  1016 
and  2.7  x  1017  cm“3,  respectively,  conditions  which  result  in 
near  Saha  values  for  ne  and  all  [Na*J. 
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SODIUM  VAPOR  DENSITY  -  [No]0,  cm'3 


Maximum  net  gain  versus  initial  sodium  density  [Na]g*  The 
numbers  in  square  brackets  are  the  electron  temperatures 
in  eV  and  the  numbers  in  parentheses  are  the  ratios  of  ne 
to  [Na]0.  F  is  the  ratio  of  ([Na*]/g*)as  to  the  Saha 
equilibrium  value. 


Absorption 


Dye  Laser 


FiRure  9  -  Schematic  of  drift  tube  and  absorption  optics  for  measurement 

of  Ne  metastable  and  resonance  atom  excitation  rate  coefficients 
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Example  of  absorption 
[Ne]  5  6,4  k  1016  cm-3 


signal  for  Ne(3P2) 
,  E/N  =  3.7  x  10" 
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Figure  11  -  Fluorescence  profile  for  X  **  588  nm  line  of  Ne  as  obtained  with 

tunable  dye  laser.  The  dashed  curve  is  the  normalized  theoretical 
Doppler  profile  including  the  isotopic  components. 
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